Wormholes in de Sitter brånes 
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In this work, we present a class of geometries which describes wormholes in a Randall-Sundrum 
bråne model, focusing on de Sitter backgrounds. Maximal extensions of the solutions are constructed 
and their causal structures are discussed. A perturbative analysis is developed, where matter and 
gravitational perturbations are studied. Analytical results for the quasinormal spectra are obtained 
and an extensive numerical survey is conducted. Our results indicate that the wormhole geometries 
presented are stable. 
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I. INTRODUCTION 

Wormholes are compact space-times with nontrivially 
topological interiors and topologically simple boundaries. 
They can be seen as connections between different uni- 
verses or topological handles between distant parts of the 
same universe. Although they are certainly exotic struc- 
tures, they appear as exact solutions of Einstein equa- 
tions with physically relevant scenarios and are compati- 
ble with the usual local physics [l[ . Samples of the work 
developed include solutions in usual general relativity [l|- 
in Gauss-Bonnet theory [(| Hl, in Brans-Dicke theory 
fl2j . and bråne world context [13l - [l6j . 
One motivation in the treatment of wormhole physics 
was due to the results of Morris, Thorne, and Yurtsever 
@; 13 j which connected time machines and traversable 
wormholes. More recently, new cosmological observa- 
tions and theoretical proposals have motivated a renewed 
interest in geometries which describe Lorentzian worm- 
holes. One of their general characteristics is that worm- 
holes must be supported by "exotic matter," which vio- 
lates usual energy conditions. Nevertheless, recent obser- 
vations suggest that the Universe may be dominated by 
some form of exotic matter (l7l . [l8| , which make worm- 
hole scenarios more plausible. Other source of geometries 
with nontrivial topology are bråne worlds. In this con- 
text, the wormhole is supported by the infiuence of a bulk 
in the bråne which describes our Universe. It is in this 
framework that the present paper is inserted. 

In this work we derive a family of asymptotically de 
Sitter wormhole solutions in a bråne world context. more 
specifically in a Randall-Sundrum-type model [19| . We 
used the effective gravitational held equations derived by 
Shiromizu, Maeda, and Sasaki [20j. As there are few sat- 
isfactory bulk solutions for compact objects in Randall- 
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Sundrum scenarios, one alternative is to build geome- 
tries in the bråne and invoke Campbell-Magaard theo- 
rems [2l[ , which guarantees their extensions through the 
bulk (locally at least). This approach has been used by 
several authors (l3l - [l6l [22J, and we will be following it 
in the present work. Moreover, global regularity of the 
bråne is expected to facilitate the construction of a reg- 
ular bulk solution [l4| . This issue will be explored in the 
space-times constructed here. 

Contrary to what has been suggested in the literature 
[Tol |23|. we obtain de Sitter solutions which are regular 
everywhere. The class of geometries studied here comple- 
ment the asymptotically flat space-times treated in (l3l — 
llq. |22 |: and the asymptotically anti-de Sitter metrics in 
[23h25| . We should mention that solutions of the effective 
Einstein equations with positive cosmological constant in 
a bråne setting were previously considered in [l6[ . While 
there is some overlap between the present work and [l6[ , 
we have explored some global issues not considered in the 
mentioned paper, such as regularity and the existence of 
cosmological horizons. As will be discussed, these points 
are particularly important for de Sitter geometries. 

If one considers the possibility of the existence of 
wormholes seriously, characteristics such as stability and 
response of wormholes to external perturbations should 
be investigated. Perturbative dynamics around worm- 
holes [26j have not been as thoroughly explored as the 
black hole problem. We further advance the perturbative 
treatment of wormhole geometries in the present work. 
Matter and gravitational perturbations are considered in 
the background of the de Sitter wormhole geometries de- 
rived here. 

The structure of this paper is presented in the fol- 
lowing. In Sec. II we have derived a family of ana- 
lytic asymptotically de Sitter solutions in an Randall- 
Sundrum-type bråne. In Sec. III the maximal extensions 
of the solutions are considered and the wormhole geome- 
tries discussed. The near extreme limit of the worm- 
hole solutions are considered in Sec. IV. Section V deals 
with the perturbative analysis of the backgrounds de- 
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rived. And finally in Sec. VI some final remarks are 
made. In this work we have used the metric signature 
diag( — h ++) and the geometric units G 4 £> = c = 1, 
where G^d is the effective four-dimensional gravitational 
constant. 

II. DE SITTER BRÅNE SOLUTIONS 



set of solutions V Cl = {(A c , B c ) , < C < d}, labeled 
by a real parameter C, such that T>d C S. 

Since Eq.(J3J) is linear in terms of B, a linear combina- 
tion of solutions with A fixed is still a solution. Moreover, 
since we are interested in deformations of the usual vac- 
uum solutions, we assume the Ansatz 

A(r)=Mr), (6) 



The basic bråne world set up is a four-dimensional 
bråne, our universe, immersed in a larger manifold, the 
bulk. It is generally postulated that the usual matter 
fields are confined in the bråne |27j . Following the ap - 
proach suggested by Shiromizu, Maeda and Sasaki [20| , 
the effective four-dimensional gravitational held equa- 
tions in the a vacuum Randall-Sundrum bråne is 



(i) 



In this effective Einstein equation, A40 is the four- 
dimensional bråne cosmological constant and E^ v is pro- 
portional to the (traceless) projection on the bråne of the 
five-dimensional Weyl tensor. Eqs. (Q} reduce to usual 
four-dimensional vacuum Einstein equations in the low- 
energy limit. 

If we impose staticity and spherical symmetry in the 
bråne, that is, 



ds 2 = -A(r)dt 2 



dr 2 
W) 

the trace of Eq. ([!]) will be 



+ r 2 (d9 2 +sm 2 0d(j) 2 ) , (2) 



i? - 4A 4 i5 , 
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where R denotes the four-dimensional Ricci scalar. The 
Eq. (J3j may be written as a constraint between the func- 
tions A and B 
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with prime (') denoting differentiation with respect to r. 

We propose to construct asymptotically de Sitter 
space-times, such that A^d > 0. In addition, we assume 
that they are "close" to the usual spherically symmetric 
(electro) vacuum solution given by A = A and B = B 0: 
with 



A (r)=B (r) = 1 



2M Q 2 A iD 2 



(5) 



where M and Q 2 are positive constants. Denoting a par- 
ticular solution by the pair (A, B) of functions which sat- 
isfy the constraint Q, we are searching for a family of 
solutions S such that: 

(i) the vacuum solution (A , B ) is an element of S: 

(ii) a generic solution (Ac 1 , Bc 1 ) G S is a continuous 
dcformation of (Aq,Bq), that is, there is (at least) one 



B(r) = B (r) - C B Hn (r) 



(7) 



with dBi in /dC = 0. Using Eqs. ([6]) and 0, the con- 
straint Q can be rewritten as a linear first order ordinary 
differential equation on the Bu n 

h(r) dBi^r)_ + ^ 
dr 

where the functions h and / are given by 

h(r) = 4A + rA' , (9) 



m = ^+*A> + 2rA»- r -^l. (10) 
r j\q 



The zero structure of h and A$ will be of great impor- 
tance. If M > 0, Q ^ and < A4D < A ext , where 
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is the critical value of A4/3, the function Aq has four 
real zeros r c , r + , r_ and r n such that r n < < r_ < 
r + < r c . Also in this region of the parameter space, 

the function h has four real zeros ro, ro-, ro and ro n 

with ron < < ro < ro- < rg. Explicit expressions 

for the several roots introduced are straightforward but 
cumbersome. Of fundamental importance in this work is 
the relation r + < ro < r c , which is always satisfied for 
< A 4D < A ext . 

The solution of (j8]) for the correction Bu n , general up 
to a multiplicative integration constant, is given by 



B lm (r)=A (r) 



(r - r — ) c 



(r-ro) Co (r-r -) Co -(r-ro„) Co " ' 

(12) 

where the positive constants cq, cq-, cq and co n are 
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(16) 



Therefore the complete solutions for A and B can be 
expressed as 



A(r) = A (r) 



B(r)=A (r) 



A 



4D 



(r c - r) (r - r+) (r — r_) (r - r„) , 

(17) 



1 - 



(r - r - 
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(18) 

It is apparent that the function B diverges in the limit 
r — > ro- Since r+ < r < r c is a natural candidate for 
the space-time static region and r+ < ro < r c , previous 
works in the literature [15|, [23| have suggested that reg- 
ular de Sitter solutions of ([3]) might not exist. However, 
we will show that this is not so. 

As will be discussed in the following sections, the main 
characteristics of this class of solutions are captured by 
the simpler case M = Q = 0. In this limit the coefficients 
r c and r can be easily expressed as 
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(19) 



and r n = -r , r 0n = -ro, c = c „ = 3/2, c _ = 1. 

The remaining constants r—, r , ro-, ro and co 

are null. The metric functions are given by 



A(r) = 1 - -o , 



B(r) =1-3 
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The energy density, radial and tangential pressures as- 
sociated with Eqs. ([2U|) and (j2"Tj) may be defined as 



(-E»)=8n 
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and are given by 
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(22) 

(23) 
(24) 
(25) 



These energy density and pressures are not generally 
positive-definite, and the effective stress-energy tensor 
(—Efy do not satisfy usual energy conditions. Still, in 
the context of this work, they should be viewed as effec- 
tive quantities, associated with a vacuum bråne model. 



III. WORMHOLES INSIDE COSMOLOGICAL 
HORIZONS 

Strictly speaking, the metric described by the func- 
tions A and B in Eqs.([l7) and ([IB]), or in Eqs. (|20| and 
PTjl , describes a space-time only for the values of the ra- 
dial parameter r such that A(r) > and B(r) > 0. The 
maximal extensions of these solutions will be presently 
treated. At this point, an important question to be 
treated is the range of the parameters for which the solu- 
tion given by Eqs. (fl"7f and (fT8)) describes an acceptable 
geometry. 

If C < 0, the functions A and B are positive for r + < 
r < r c . But r + < ro < r c , so B in Eq. (|18p is divergent 
at ro- The geometry is well-defined and static for r > ro, 
but its curvature invariants are not bounded, as seen by 
the behavior of the Kretschmann scalar near ro 



lim |f? Qj a 7 5i? 



a/?7<5 I 



— > OO 



(26) 



Therefore, for this case a naked curvature singularity is 
present at r — > ro- This solution will not be further 
explored in the present work. 

If C — 0, we recover the usual Reissner-Nordstrom-de 
Sitter vacuum solution, and the regular region is given by 
r + < r < r c . As is well-known, in nonextremal regimes 
the surfaces r — r + and r — r c describe an event and an 
cosmological horizons in the maximal extension, respec- 
tively. One interpretation for this result is that, although 
the solutions with C ^ and the Reissner-Nordstrom-de 
Sitter black holes have very different global characteris- 
tics, they nevertheless are locally arbitrarily close. 

If C > 0, the function B is not positive-definite be- 
tween r + and r c . It has a third zero at r = r t hr- The 
relevant point is 



r + <r Q < r t hr < r c 



(27) 



with the functions A and B positive-definite and analytic 
for r t hr < r < r c . Therefore, the chart (t,r,9,<j>) is valid 
in the region r t hr < r < r c . The analytic extension 
beyond r = r t hr is suggested with the use of the proper 
length £ as radial function, where 



d£(r) 
dr 



1 



VW) 



(28) 



Choosing an appropriate integration constant in Eq. (|28|) . 
the region r t hr < r < r c is mapped into < £ < £ m ax, 
with a finite l ma x- The extension is made analytic con- 
tinuing the metric with —i max < £ < £ m ax- The result- 
ing geometry has a wormhole structure, with a throat at 
r = r th r- 
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Figure 1. Conformal diagram of wormhole the solution inside 
a cosmological horizon. Dashed line denotes the wormhole 
throat. 



The extension beyond r — r c can be made, for ex- 
ample, with the ingoing and outgoing Eddington charts 
(u,t,9,4>) and (v,t,0,(j>), where u, v are the light-cone 
variables 



u = t — r* and v = t + r* 



(29) 



The radial variable r* is the tortoise coordinate, defined 
as 



dr 



(30) 



In the maximal extension, the surface l = £ rnax (r = r c ) 
is a Killing horizon, interpreted as a cosmological hori- 
zon. A more physical interpretation of the geometry is a 
spherically symmetric wormhole inside an exponentially 
cxpanding Universe. 

Charts based on the tortoise coordinate r* or the 
u — v coordinates can be used to cover all the static 
region. In this case, {(t, é, 9, (f) , l G (-l max , +i m ax)} 
is mapped into {(t, r*, 9, (f>) , r+ G (— co, +oo)} or 
{(u, v,9,<fi) ,u E (— oo, +cxd), v G (— oo, +oo)}. These 
coordinate systems will be used in the perturbative 
analysis of the wormhole. 

Applying standard procedures (see, for example, [28]), 
the Penrose diagram of wormhole geometry can be ob- 
tained. This diagram is present in FigfT] 



IV. THE NEAR EXTREME LIMIT 

One limit where the geometry is simpler and its per- 
turbative analysis is much easier (as we will see in the 
next section) is the near extreme regime. We will treat 
this limit in the present section. 

In order that the background characterized by the met- 
ric functions (|17l) and (|18[) describes a wormhole, the real 
parameter C cannot be arbitrarily large. As C grows, 
rthr approaches r c . The extreme value for C (C ex t), such 



Tthr 



r c with < C < C ex t , is given by 



C' 



(Xc - r ) c ° (r c - r _) Co (r c - r 0n ) c 
(r c -ro-) co - 



(31) 



We will consider in this section the near extreme limit 
case, where C < C ex t, that is, C is very close (but still 
smaller) to the maximum value C ex t- So, it is natural to 
define the dimensionless parameter 



rth, 



ro 



(32) 



With this defmition, < S < 1, since r < r t h r < r c . 
The near extreme regime can be characterized in terms 
of 5 as the limit < 5 <C 1. In fact, it can be shown that 
C/C ext = l-0(5). 

In the near extreme limit, the metric functions A and B 
can approximated by the linear and quadratic functions 
A n ~ ext and B n ~ ext respectively, 



A(r) 



Hr) 



(r c - r) 



(33) 



B(r) ra B n ~ ext (r) = B^ ext (r - r thr ) (r c - r) , (34) 



with the positive constants Aq 

An-ext _ -A-4.D 

Ar 



and B™~ ext 



given by 

2 (rc - r+) (r c - r_) (r c - r„) , (35) 
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(36) 



It is important to stress that the causal structure of the 
space-time is not modified in the near extreme limit. The 
geometry still describes a wormhole inside cosmological 
horizons, with its Penrose diagram shown in Fig[TJ One 
important geometrical quantity is the surface gravity at 
the cosmological horizon. In the near extreme limit it 
can be explicitly calculated in terms of the roots of A 
and h 



d^A{r)B{r) 



dr 



a n—ext r>n — ext 

A o B o 



(r c - r ) 5 1 ' 2 . (37) 



Although the cosmological horizons may be seen as 
close in the near extreme limit, this is not necessarily so. 
In fact, the proper radial distance between the horizons 
can be arbitrarily large even in the near extreme regime. 
Taking the case M = Q = for simplicity, the max- 
imum value for the radial proper distance (£ m ax), half 
the proper distance between the two cosmological hori- 
zons, is lmax = iTir c /2, which can be arbitrarily large 
as r c — > oo — > 0). The proper distance between 
horizons is then unbounded. 
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As discussed in previous and following sections, charts 
based on the tortoise coordinate introduced in Eq. tptO)) 
are very convenient for several applications. In the near 
extreme regime, the metric can be explicitly written in 
terms of r* 



ds 2 = 6A£- ext (r c - r ) sech 2 { Kc r+) (-dt 2 + dr 2 ) 

+ [r c — S (r c — ro) sech 2 (fc c r*)] 2 d£l\ . (38) 



V. PERTURBATIVE DYNAMICS AND 
STABILITY ANALYSIS 

A. General considerations for the perturbative 
treatment 

Once the background geometry is established, one next 
step is to determine its response under small perturba- 
tions. In the lowest order, background reaction can be 
ignored, and the dynamics is restricted to the matter 
and gravitational perturbations in a fixed geometry. As 
a prototype of matter, we will consider massless and mas- 
sive scalar held, not necessarily minimally coupled to the 
background. The gravitational perturbation analysis will 
be limited here to the axial mode dynamics. 

A massless scalar perturbation held <f> is characterized 
by the Klein-Gordon equation 



□$ = 



(39) 



Decomposing the scalar held $ in terms of an expansion 
in spherical harmonic components 



(40) 



Lm 



the Klein-Gordon equation give us a set of decoupled 
equations in the form 



d 2 ipi d 2 i\)i 
dt 2 dr 2 



V sc (r(r*))iPi , 



(41) 



labeled by the multipole index l, with l = 0, 1, 2, The 

tortoise coordinate r* was introduced in Eq. (|30[) , and Yi m 
denotes the spherical harmonic functions. Using results 
in [29] , the scalar effective potential is expressed in terms 
of r as 

V sc (r) = l -^-^A + -A A' Q - £ (A B Hn )' . (42) 

Typical profiles for the scalar effective potential are pre- 
sented in Figf2] 



We will consider gravitational perturbations in the 
bråne geometry following the treatment in (29| . In gen- 
eral, the gravitational perturbations depend on the tidal 
perturbations, namely, first-order perturbations in 





-50 -25 



-50 -25 25 50 



Figure 2. Scalar and gravitational effective potentials (right 
and left panels, respectively) in terms of the tortoise coor- 
dinate r*. The wormhole parameters used in the plots were 
A 4D = 0.01, M = 1.0, Q = 0.5 and S = 0.7 (r thr = 14.20 and 
r c = 16.23). 



{8E^ V ). Since the complete bulk solution is not known, 
we shall use the simplifying assumption 8E^ V = 0. This 
assumption can be justihed at least in a regime where 
the energy carried in the perturbation processes does not 
exceed the threshold of the Kaluza-Klein massive modes 
[27| . Analysis of gravitational shortcuts [3(1 [H| also sup- 
ports this simplification, suggesting that gravitational 
fields do not travel deep into the bulk. Within these 
premises, the gravitational perturbation equation is 



SR^ = 



(43) 



Following [23|, the gravitational axial perturbations 
are given wave functions Zi , satisfying a set of equations 
of motion with the form (|45|) , labeled by a multipole in- 
dex l (l = 2, 3, 
given by 



.). The effective potential in this case is 



Vgrav (r) 



-C 



g-i) , , 2 2 

~2 A ° + ^2 A - 

"2 \ 2 

— A Bi in — — (AoBlin) 



-A A 

r 



(44) 



Typical profiles are presented in Fig(5] 

Of particular interest in the perturbative dynamics are 
the so-called quasinormal mode spectra. Consider a wave 
function R, in the present case the scalar or gravitational 
perturbation (■ø; or Zi), subjected to an effective poten- 
tial V (V sc or Vgrav)- The quasinormal modes are Solu- 
tions of the "time-independent" version of Eq. (|4Tj) , 



d 2 R. 
drl 



+ (lj 2 -V)R lj = 



(45) 



satisfying both ingoing and outgoing boundary condi- 
tions asymptotically: 



lim R u e 

r t — f =poo 



= 1 



(46) 
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The "frequency domain" wave function R^ associated 
with a givenquasinormal mode lo is given by the Laplace 
transform 1 321 of the function R as 



Rit^e^dt 



(47) 



with lu extended to the complex plane. 

The relevance of the quasinormal mode calculation is 
twofold. They determine the dynamical evolution of the 
wave function when the wave equation is subjected to 
bounded initial conditions. Moreover, lm (to) < is a 
necessary condition for the stability of the perturbation. 
The determination of quasinormal mode spectra for the 
wormhole geometry will be made with analytical and nu- 
merical techniques in the next subsections. 

As will be discussed in the following, the dynamics of 
the perturbations considered can be analytically treated 
in the near extreme limit introduced in Sec. IV. Beyond 
this regime, numerical tools are necessary. In order to 
analyze quasinormal mode phase and late-time behavior 
of the perturbations, we apply a numerical characteristic 
integration scheme based in the light-cone variables u and 
v in Eq. ([29|) . used, for example, in (33l - l36| . 



B. Spherically symmetric scalar mode (l = 0) 

The scalar held perturbation has a spherically sym- 
metric (l = 0) mode. This mode is distinct because its 
associated effective potential is not positive-definite, as 
illustrated in FigJ^l This point raises the question of 
whether the time evolution of the scalar held is stable. 
One important result of this work is that, in our exten- 
sive numerical investigation, the perturbation is always 
bounded, that is, no unstable modes were observed. 

The presence of relevant negative peaks in the scalar 
potential with l = is a potential complication for the 
calculation of the quasinormal frequencies. Nevertheless, 
the direct integration scheme used in (33l - 436j can be suc- 
cessfully employed in the present case. We will discuss 
in the following some important points observed in the 
scalar held evolution in the wormhole background con- 
sidered. 




Figure 3. Scalar field evolution with l — and r* = 0, for 
several values of S. For the wormhole geometries considered, 
the parameters A^d = 0.01, M = 1.0 and Q = 0.5 (r c = 
16.23) were used. 
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Figure 4. Dependence of the asymptotic value for the l = 

scalar mode (Vo°') with the parameter S. The hullets in- 

dicate numerical results, and the dashed line denotes a å 2 
power-law. For the wormhole geometries considered, the pa- 
rameters A 4 d = 0.01, M = 1.0 and Q = 0.5 (r c = 16.23) were 
used. 



A nonusual feature observed in the scalar dynamics is 
that the field ip, for a fixed value of r*, tends to a non- null 
constant '00°'' f° r large t: 

lim ^ =0 -► 4 0) • (48) 

t— ¥00 

This point is illustrated in FigJ3] A similar qualitative 
behavior was observed in other de Sitter geometries (35l . 
& 

In the near extreme regime, the late-time field evo- 
lution can be better explored. The intermediate- and 
late-time field evolution has the form 

tø=o^° ) +^ 1) e-"' t > (49) 



with 

4 0) oc S 2 . (50) 

and k c denoting the surface gravity at the cosmological 
horizon, calculated at Eq. (|37p in the near extreme limit. 
The dependence of -øq -* with the parameter å is illus- 
trated in FigH] 



7 



10 
10 * 

: lo" 12 
lo 16 
io' 20 



5 = 0.6 
5 = 0.8 
5 = 0.95 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 5. Scalar field evolution with l — 1 and r* = 0, for 
several values of S. For the wormhole geometries considered, 
the parameters A4D = 0.01, M = 0.1 and Q — 0.05 (r c — 
17.22) were used. 



C. Higher multipole modes (l > 0) 

A general feature of the effective potentials considered 
when l > is that they are positive-definite. This point 
implies that the dynamics is always stable for non-null 
l. Other relevant characteristics of both potentials are 
the typically complicated profiles near r* w 0, as illus- 
trated in Fig[2j This latter point makes the WKB-based 
methods in [37-39] not effective in the present 
explicitly checked by us. The direct integration schemes 
used in [33h36] can still be successfully employed. Ana- 
lytic results will be available in the near extreme regime. 



For the scalar perturbation with 1 = 1, the main qual- 
itative characteristics of its perturbative dynamics are 
described as follows. If å is close to 1 (C/C ex t small), the 
late-time decay is (nonoscillatory) exponential, 



i>l=i 



(51) 



with k < 0. This result is consistent with the scalar dy- 
namics around other asymptotically de Sitter geometries 
[H, Ull- We illustrate this result in Fig[5] For smaller 
values of 5 (larger C/C ex t), the decay is oscillatory, with 
an exponential envelope, 



7=1 



Im(wQ C ) t -iKe(uig C ) t 



(52) 



where u!q c is the fundamental (lowest absolute value of 
the imaginary part) quasinormal frequency associated 
with the l = 1 scalar mode and lm (uJq c ) < 0. We il- 
lustrate this result in Fig. [5] 

Typical profiles for the dependence of the parameters 
fe, lm (uJq c ) and Re (wq c ) on 6 are shown in FigJBl From 
these results, we see that the shift of oscillatory and 
nonoscillatory modes at t — > 00 is determined by the 
relative magnitude of k and Im(wQ C ). If |fe| > |Im(wQ C )| 
(small S) , the nonoscillatory mode is suppressed for large 
t, and the oscillatory phase dominates. If |Im (wq c )| > \k\ 
(large enough S), the oscillatory mode is suppressed, and 
a late-time nonoscillatory decay is observed. 
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Figure 6. Dependence of the exponential coefficient k and the 
scalar fundamental quasinormal frequency u;o c on the param- 
eter S, for the 1 = 1 mode. For the wormhole geometries con- 
sidered, the parameters A4D = 0.01, M = 0.1 and Q = 0.05 
(r c = 17.22) were used. 



For scalar or gravitational perturbations with l > 1 the 
intermediate and late-time dynamics is dominated by an 
oscillatory exponential decay. The scalar and gravita- 
tional perturbations can be well characterized by their 
fundamental quasinormal frequencies (wq c and ujQ Tav ): 



^~e- iw °*, (53) 

Z l ~ e -^o rav t . (54) 

These results are illustrated in FigJTl We have not ob- 
served nonoscillatory exponential decays for the scalar or 
gravitational perturbations with l > 1, considering values 
of C/C ext as low as 1CT 4 (8 w 0.999). 



In the near extreme regime (0 < å <C 1 or C/C ex t < 1), 
considered in Sec. IV, the scalar and gravitational quasi- 
normal mode spectra can analytically determined. Ex- 
plicit analytic expressions for the functions V r sc (r(r*)) 
and V r sc (r(r Jr )) are usually not available, except in par- 
ticular limits. One of these limits is the near extreme 
regime. Following an approach similar to the one used in 
[40,|4l|, the result (f3"5)) allows both effective potentials to 
be written as 

V(r(r*))= ^ ax (55) 
cosh («er*) 

with the surface gravity k c presented in Eq. (|37j) . The 
constants V m ax, for the scalar and gravitational pertur- 
bations (V^ax an d Vmax ■> respectively) are 



ysc = 

' max 



5A iD l{l + l) 



(r c - r ) (r c - r+) (r c - r_) (r c - r n ) 



3ri 
with l > , 



(56) 



s 



Table I. Fundamental quasinormal frequencies for the scalar perturbation for several values of 8 and l. For the wormholc 
geometries considered, the parameters A4/5 = 3 and M = Q = (r c = 1) were used. Relative errors (A%) for the near extreme 
results are indicated. 







Direct intej 


jration 


Near extreme results 


1 




Re (wg c ) 


lm (uj a c ) 


Re(wo c ) (A%) 


lm (wg )(A%) 


1 


0.001 


0.01660 


-0.02143 


0.01659 (0.06%) 


-0.02142 (0.05%) 


1 


0.01 


0.05224 


-0.06806 


0.05246 (0.42%) 


-0.06773 (0.49%) 


1 


0.1 


0.1580 


-0.2175 


0.1659 (5.00%) 


-0.2142 (1.52%) 


1 


0.3 


0.2318 


-0.4130 


0.2874 (24.0%) 


-0.3710 (10.2%) 


1 


0.5 


0.2084 


-0.5747 


0.3710 (78.02%) 


-0.4789 (16.7%) 


2 


0.001 


0.04177 


-0.02143 


0.04175 (0.05%) 


-0.02142 (0.05%) 


2 


0.01 


0.1321 


-0.06806 


0.1320 (0.08%) 


-0.06773 (0.48%) 


2 


0.1 


0.4192 


-0.2248 


0.4177 (0.36%) 


-0.2143 (4.67%) 


2 


0.3 


0.7256 


-0.4307 


0.7232 (0.33%) 


-0.3710 (13.8%) 


2 


0.5 


0.9157 


-0.6025 


0.9336 (1.95%) 


-0.4789 (20.5%) 


2 


0.7 


1.017 


-0.7086 


1.105 (8.65%) 


-0.5667 (20.0%) 


2 


0.9 


1.012 


-0.9078 


1.253 (23.8%) 


-0.6426 (29.2%) 



Table II. Fundamental quasinormal frequencies for the gravitational perturbation for several values of 8 and l. For the wormhole 
geometries considered, the parameters A413 = 3 and M = Q = (r c = 1) were used. Relative errors (A%) for the near extreme 
results are indicated. 



Direct integration Near extreme results 



l 


S 


ReK™") 


ImK™ ) 


Re (w*™*) (A%) 


ImK 9m ") (A%) 


2 


0.001 


0.03178 


-0.02141 


0.03177 (0.03%) 


-0.02142 (0.05%) 


2 


0.01 


0.1004 


-0.06758 


0.1005 (0.09%) 


-0.06773 (0.22%) 


2 


0.1 


0.3167 


-0.2095 


0.3177 (0.32%) 


-0.2142 (2.24%) 


2 


0.3 


0.5457 


-0.3354 


0.5503 (0.84%) 


-0.3710 (10.6%) 


2 


0.5 


0.6962 


-0.4383 


0.7104 (2.04%) 


-0.4789 (9.26%) 


2 


0.7 


0.8054 


-0.5074 


0.8405 (4.36%) 


-0.5667 (11.7%) 


2 


0.9 


0.8647 


-0.5622 


0.9531 (10.2%) 


-0.6426 (14.7%) 


3 


0.001 


0.05669 


-0.02142 


0.05668 (0.02%) 


-0.02142 (0.00%) 


3 


0.01 


0.1792 


-0.06778 


0.1792 (0.0%) 


-0.06773 (0.07%) 


3 


0.1 


0.5664 


-0.2158 


0.5667 (0.05%) 


-0.2142 (0.74%) 


3 


0.3 


0.9785 


-0.3803 


0.9815 (0.31%) 


-0.3710 (2.24%) 


3 


0.5 


1.254 


-0.5012 


1.267 (1.04%) 


-0.4789 (4.45%) 


3 


0.7 


1.458 


-0.6074 


1.499 (2.81%) 


-0.5667 (6.70%) 


3 


0.9 


1.577 


-0.7013 


1.700 (7.80%) 


-0.6426 (8.37%) 



yg 



rav 
max 



8A iD (J + 2)(J-1) 



(r c - r ) (r c - r+) (r c - r_) (r c - r n ) 



3rj 



with / > 1 



(57) 



The potential in (|55|) is the so-called Poschl- Teller po- 
tential [42]. It has been extensively studied, and, in par- 
ticular, the quasinormal modes associated with it have 
been calculated [H, Bj. Using the results in [431. E3|. we 
have for the scalar and gravitational quasinormal mode 
spectra, in the near extreme regime 



iysc 1 
r max 



(58) 



(59) 



with k c 



and VZ a a x v 



given by expressions (|3T)) . 



V s 

■ m 

and (|57|) respectively. 

The fundamental {11 = 0) modes dominate the late- 
time decay. We stress the excellent concordance of the 
analytical expressions (|55)l - ((55|) with the numerical re- 
sults in the near extreme regime. Moving away from the 
near extreme limit, we consider the quasinormal spectra 
for higher values of 8. A direct integration approach has 
been used. Quasinormal frequencies for the scalar and 
gravitational perturbations are presented in Tables|T]and 

m 

In all numerical calculations performed, the concor- 
dance between the numerical and near extreme approxi- 
mation improves as 8 is made smaller. This is a consis- 
tency check for the numerical results and an indication 
that the near extreme results are indeed adequate when 
the appropriate limit is tåken. We illustrate this point 
in Tables Hl and UTI Moreover, the analytical expression in 
Eq. (jnU) for the gravitational sector appears to work well 
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Figure 7. Scalar and gravitational perturbations with / = 2 
and = 0, for several values of S. For the wormhole geome- 
tries considered, the parameters A4D ~ 3 and M = Q = 
(r c = 1) were used. 



even when the condition 5 <C 1 is not strictly satisfied, as 
suggested by the data presented, for example, in Table 



VI. FINAL REMARKS 

We have obtained a family of exact solutions of the 
effective Einstein equations in a asymptotically de Sit- 
ter Randall-Sundrum bråne. This family includes naked 
singularities, but also solutions which describe worm- 
holes. Maximal extensions of the solutions were studied. 
We have shown that the extensions describe Lorentzian, 
traversable, wormhole space-times which connect regions 
bounded by cosmological horizons. It should be noted 
that, although the existence of a local, asymptotically 
de Sitter, solution for a metric in a Randall-Sumdrum 
scenario might be expected, it is not obvious that there 



would exist solutions regular everywhere. The explicit 
solutions constructed here have this characteristics. 

One basic requirement, if the geometries obtained are 
to be considered as physically relevant, is the stability 
of the derived geometries under first-order perturbations. 
We have treated this question here considering scalar and 
axial gravitational perturbations. An important result in 
the perturbative analysis performed in this work is that 
no unstable modes were found. 

Moreover, the detailed numerical and analytical treat- 
ment presented sketches a picture of the perturbative 
dynamics. Scalar spherically symmetric modes typically 
decay to a nonzero constant asymptotically. This is rem- 
iniscent of a feature already observed in considerations 
involving de Sitter black holes [34l . [35[ . Although oscilla- 
tory and nonoscillatory decays bounded by exponential 
envelopes were observed, no power-law tails appeared, 
which also resembles the dynamics around asymptoti- 
cally de Sitter black holes [34 , [35( . 

One interesting limit of the geometries derived in this 
work is their near extreme regime. This limit is inter- 
esting because the geometry becomes very simple, while 
still preserving the causal structure of the nonextreme 
case. In fact, in the near extreme regime the quasinor- 
mal spectra of the perturbations considered can be an- 
alytically determined, which is something not common 
in the literature. Moreover, the comparison between the 
full numerical results and the near extreme approxima- 
tion shows good agreement for the fundamental overtone. 
We consider this result a strong argument for the val- 
idation of both approaches. Besides, the near extreme 
analytical results appear to describe reasonably well the 
gravitational quasinormal spectra even outside this limit. 
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